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We explore an alternative mechanism for the production of gravitino dark matter whereby relic 
gravitinos originate from the decays of superpartners which are still in thermal equilibrium, i.e. via 
freeze-in. Contributions to the gravitino abundance from freeze-in can easily dominate over those 
from thermal scattering over a broad range of parameter space, e.g. when the scalar superpartners 
are heavy. Because the relic abundance from freeze-in is independent of the reheating tempera- 
ture after inflation, collider measurements may be used to unambiguously reconstruct the freeze-in 
origin of gravitinos. In particular, if gravitino freeze-in indeed accounts for the present day dark 
matter abundance, then the lifetime of the next-to-lightest superpartner is uniquely fixed by the 
superpartner spectrum. 



I. INTRODUCTION 

Supersymmetry is an elegant and well-motivated ex- 
tension of the standard model which solves the hierarchy 
problem and carries extensive phenomenological conse- 
quences. Despite its successes, however, supersymme- 
try suffers from an assortment of cosmological difficulties 
which are referred to collectively as the "cosmological 
gravitino problem" , which has two components: 

a) Late decaying superpartners can produce electro- 
magnetic or hadronic radiation that can adversely 
affect Big Bang Nucleosynthesis (BBN). 

b) Relic gravitinos produced by scattering and decay- 
ing superpartners can overclose the universe. 

The gravitino problem highlights a tension between su- 
persymmetry and cosmology which is highly robust. This 
is so because the existence of the gravitino is required by 
supergravity, and because the couplings of the gravitino 
to superpartners are uniquely fixed by soft masses and a 
single additional parameter, the gravitino mass. 

At the same time, there are a number of approaches 
by which to address these issues. For example, a) may 
be evaded if the scale of supersymmetry breaking is very 
high, as in anomaly mediation, in which case the grav- 
itino will decay safely before BBN. Alternatively, super- 
partner decays to the gravitino can be made sufficiently 
rapid if the scale of supersymmetry breaking is low or 
intermediate, as in gauge mediation. 

Likewise, b) can be resolved if ^ keV, in which 

case the gravitino is simply too light to overclose the 
universe^. For 1113/2 > keV, the authors of [l| famously 
showed that b) can be avoided if the reheating temper- 
ature after inflation, Tr, is below a critical value which 
depends on the superpartner spectrum and is shown in 
Fig. [1] Numerous new physics proposals — for example. 



ones including new stable charged particles or super- 
weakly interacting particles — suffer from an analogous 
overclosure problem which may be evaded by appropri- 
ately lowering Tr. 

The conventional wisdom is that low Tn is disfavored; 
for example, > 10^ GeV for high scale leptogenesis Q. 
As a consequence, the vast majority of papers on grav- 
itino cosmology have focused on the portion of Fig. [T] at 
high Tfj. However, the cosmological baryon asymmetry 
can be generated at much lower temperatures, for exam- 
ple via soft leptogenesis and in this case gravitinos 
become a virtue rather than a problem: not only are 
both a) and b) resolved, but gravitinos can fully account 
for the observed dark matter for a wide range of masses, 
as depicted in Fig. [TJ 

In this paper, we investigate a more or less overlooked 
regime of gravitino cosmology corresponding to the ver- 
tical incline in critical shown as a function of 
in Fig. [1] Here gravitino dark matter arises dominantly 
from the freeze-in mechanism, which was studied in some 
generality in In this setup, a feebly interacting dark 
matter particle is produced via the decays of particles 
which are still in thermal equilibrium. Crucially, since 
the decay rates of these particles fix the final abundance 
of dark matter, the associated lifetimes are hence con- 
strained by the observed dark matter abundance. Fur- 
thermore, because the production is dominated at low 
temperatures, the freeze- in abundance is largely indepen- 
dent of Tr, explaining the vertical incline in critical T^. 
For the case of gravitino freeze-in, decays of superpart- 
ners in thermal equilibrium produce a final yield of 
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where i sums over all superpartners, rui and Fj are super- 
partner masses and partial decay widths to the gravitino, 
and mpi is the reduced Planck mass^. 



^ However, note that this class of theories is in tension with warm 
dark matter constraints. 



^ Throughout, sums over superpartners will implicitly include a 
degeneracy factor — for instance, a factor of 8 for gluinos, etc. 
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FIG. 1: Contours of ^3/3/1^ = 0.11 for gaugino masses fixed 
to {mi, ma, nig} = {100,210,638} GeV. The {red, orange, 
yellow, green, blue} contours correspond to universal scalar 
masses {500 GeV, 1 TeV, 2 TeV, 4 TeV, 8 TeV}. 

As we will see, the freeze-in abundance of gravitino 
dark matter depends solely on the superpartner spectrum 
and ?Ti3/2, a quantity which is straightforwardly inferred 
from the mass and lifetime of the NLSP when it decays to 
the gravitino LSP. Thus, for a given superpartner spec- 
trum, the constraint of fty'^^^h'^ ~ 0.11 entirely fixes the 
lifetime of the NLSP. Because these quantities are exper- 
imentally accessible, we chance upon the rather amazing 
prospect of reconstructing the origin of gravitino dark 
matter through collider measurements. For example, for 
degenerate heavy squarks and ghiinos at a mass m, the 
NLSP lifetime is 

in-7 f "^NLSP \ f m Y fn\ 

TNLSP - 10 sec I , 2) 

V300 GeV/ \mNLSP / 

if gravitino freeze-in accounts for the present day abun- 
dance of dark matter. Note that this proposal is a spe- 
cific instance of the generalized cosmological scenario dis- 
cussed in 0, . 

While the lifetime tnlsp indicated by Eq. ([2]) is ef- 
fectively long-lived on collider time scales, a number 
of theoretical and experimental collaborations have sug- 
gested that the LHC is capable of measuring the long- 
lived decays of the sizable number of charged or colored 
metastable NLSPs which will typically slow and eventu- 
ally stop within the detector material. Sufficiently long 
lifetimes can easily arise in theories of split supersym- 
metry fj^, as well as theories with very weakly coupled 
particles like gravitinos Q, axinos goldstini [l^, □J, 
sterile sneutrinos [l3|, and dark matter @, [l^. Hence, 
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FIG. 2: Contours of Q^/2h^ = 0.11 for universal scalar masses 
fixed to 500 GeV. The {red, orange, yellow, green, blue} con- 
tours correspond to a bino mass mj; — {500 GeV, 1 TeV, 2 
TeV, 4 TeV, 8 TeV}, where and rrig are fixed assuming 
gaugino mass unification at Mqvt ~ 10^^ GeV. 

stopped NLSPs allow for a range of 10~® — 10^ sec to 
be probed in early LHC running, and indeed bounds on 
stopped gluinos have already been set by the CMS col- 
laboration [l3l- At higher luminosities, neutral NLSPs 
might also be probed if their lifetimes lie in the range 
10~^ — 10~^ sec. As such, gravitino freeze-in offers a 
novel mechanism of dark matter generation which has 
direct implications for the LHC in the near term. 



II. GRAVITINO COSMOLOGY 

Assuming that the messenger scale of supersymme- 
try breaking is below the Planck scale, then the grav- 
itino is the lightest of all the superpartners and is thus 
an attractive R-parity stabilized dark matter candidate. 
Typically, the gravitino mass is considered in the range 
keV < m3/2 ^ 1 GeV, where the lower bound arises 
from warm dark matter constraints and the upper bound 
arises from tension with BBN^. Broadly speaking, grav- 
itinos are produced via three distinct physical mecha- 
nisms, each with a much different dependence on the re- 
heating temperature after inflation, Tr, and the gravitino 
mass, 7713/2. 



^ The quantitative BBN bound on varies with the nature 

and mass of the NLSP. Moreover, in some cases it can be evaded 
altogether, e.g. with sneutrino NLSP or R-parity violation. 
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FIG. 3: Prediction of tklsp for gravitino dark matter arising 
from freeze-in, for rriNLSP = 300 GeV. The sum over super- 
partner contributions is assumed to be dominated by gluinos 
and degenerate squarks as shown in Eq. ((8)|. 



A. Modes of Production 

First, there is a contribution to the gravitino abun- 
dance arising from NLSPs vi^hich freeze-out and then de- 
cay to the gravitino, as in the so-called superWIMP sce- 
nario [l^. This contribution is highly model dependent 
and can easily be negligible since the final relic gravitino 
abundance is down by a factor of m3/2/m^NLSP relative to 
the freeze-out abundance of the NLSP, which can itself 
be small if the NLSP has strong annihilation channels. In 
addition, BBN is in tension with the superWIMP mecha- 
nism as the origin of the dark matter [l5| . For these rea- 
sons, we will ignore superWIMP contributions and focus 
on other gravitino production mechanisms. 

A second source of gravitino production is the ther- 
mal scattering of superpartners in the early universe. 
Since the goldstino couples to gauginos through a dimen- 
sion five operator, the associated scattering processes are 
dominated at high temperatures, and so the final abun- 
dance of gravitinos depends linearly on Tr. For example, 
the yield of gravitinos from gaugino scattering goes para- 
metrically as 



3/2 



(3) 



where a — 1,2,3 sums over the gauge group and nia 
are the gaugino masses. This scattering contribution has 
been the primary focus of existing work on gravitino dark 
matter, and corresponds to the straight, sloped portions 
of the contours in Figs. [T] and [U which depict contours 
of the total gravitino abundance f23/2^^ = 0.11 in the 



(7713/2, Tii) plane for different choices for the superpartner 
spectra. 

Lastly, gravitinos may be produced by freeze-in: that 
is, from the decays of superpartners which are still in 
thermal equilibrium. The near vertical portions of the 
curves in Figs. [T] and [5] correspond to freeze-in, which is 
infrared dominated and thus independent of T^. Plug- 
ging in for Ti in Eq. ([1]) yields the parametric dependence 



Trdecay 
-^3/2 



(4) 



where i sums over all superpartners. As T/j drops be- 
low the superpartner masses, the superpartners are not 
efficiently produced from reheating, and hence gravitino 
production arises from the exponentially tiny Boltzmann 
tail fixing the abundance of superpartners. 

The size of the freeze-in region in Figs. [1] and [2] is de- 
termined by the competition between ^3^/2" ^^'^ ^3/2^^ ■ 
Thus, let us define to be the reheating temperature at 
which these two quantities are equal. Clearly, the range 
in Tii in which freeze-in is operative runs from the super- 
partner masses up to T^. Since this cross-over value of 
the reheating temperature goes as 
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(5) 



this implies that the freeze-in region will diminish for 
larger gaugino masses, and will grow for larger scalar 
masses. 

The trend implied by Eq. ([5]) is verified m Fig. m 
which depicts 513/2^^ = 0.11 from a total gravitino 

yield of ^3/2 - ^3^2"" + ^3%'"'' in the (7713/2, Tfl) plane. 
One can see that the freeze-in region grows with in- 
creasing scalar masses. This is the case because heav- 
ier scalars imply larger decay rates without commen- 
surately larger contributions from thermal scattering at 
high temperatures. For the allowed region of gravitino 
masses, 1 keV < 7713/2 < 1 GeV, the reheating tem- 
perature required for gravitino dark matter varies over 
100 GeV < Tfl < 10^ GeV. 

We now stress a key point: a sizable fraction of this 
range corresponds to gravitino freeze-in and is thus in- 
sensitive to Tfi, which is fortuitous because Tr is not 
an experimentally accessible quantity. Hence this allows 
for the unique possibility of reconstructing the freeze- 
in origin of gravitino dark matter from LHC measure- 
ments. As shown in Fig.dl for squarks accessible at LHC 
{rUg < 2 TeV) this occurs in about 30% - 50% of the 
logarithmic range of Tr, while for heavier squarks the 
range is even greater. From a theoretical standpoint, 
it is straightforward to make the scalars quite heavy 
while keeping the gauginos light with an R-symmetry. 
Note that while Fig. [1] was produced assuming degener- 
ate scalar masses, an almost identical plot results if the 
top and bottom squarks are pushed down to near the 
weak scale, as considered in fig . 



4 



Fig. [2] also verifies the trend indicated by Eq. ([5]), since 
it shows the freeze-in region diminishing for increasing 
gaugino masses. This is the case because heavier gaugi- 
nos imply larger scattering cross-sections at high energies 
and thus a larger contribution arising from Yy^^^. From 
a top-down viewpoint, theories with very heavy gauginos 
and light scalars are difficult to accommodate, since very 
large values of rrig tend to drag up rriq and exacerbate 
fine-tuning of electroweak symmetry breaking. Thus, we 
conclude that it is actually theoretically difficult to ob- 
tain theories in which the freeze-in region is small, and 
so a large freeze-in region is a typical feature of many 
reasonable models. 

Also, let us not that immediately after freeze-in, the 
produced gravitinos are highly relativistic, but they be- 
come non-relativistic as the temperature drops below 
their mass, yielding cold dark matter. At the same time, 
very light gravitinos, 7713/2 < 10 keV, yield warm dark 
matter, especially since freeze-in arises from decays of 
superpartners in the exponential tail of their thermal dis- 
tribution. 



B. Reconstructing the Origins of Dark Matter 

Let us now consider the extent to which the freeze-in 
origin of gravitino dark matter might actually be verified 
at the LHC. Assuming that the present day abundance 
of gravitinos arises entirely from freeze-in, 1^3/2 = ^3/2^^^ ■• 
one can rewrite Eq. ([T|) as 



TO3/2-!^3/2 



Y.J. = 



0.26 

3/2 



which is obtained from 
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while normalizing the partial widths of the superpartners 
decaying into gravitinos, F^, with respect to the NLSP 
decay width, Fnlsp, so Ti/mf = FNLSp/m-NLSP- From 
Eq. ^ we see that ^13/2^^3/2 oc is dominated by 

the very heaviest superpartners. Assuming that the su- 
perpartner spectrum is measured, Eq. ([6]) can be inverted 
to yield a critical prediction for the NLSP lifetime in 
terms of the superpartner spectrum: 
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The second line corresponds to an approximation in 
which the gravitinos are produced dominantly by squarks 



and gluinos. Also, while the value of (/* actually varies 
substantially with temperature, for this approximate ex- 
pression we have normalized its value to 150, which lies 
somewhere between the for the standard model and 
the minimal supersymmetric standard model. The corre- 
sponding prediction obtained by numerically solving the 
Boltzmann equations is shown in Fig. |31 

Thus, TNLSP may be as short as 10^^ sec for a squashed 
supersymmetric spectrum and as long as 100 sec, the ap- 
proximate bound from BBN, for an extremely split spec- 
trum. Note that this entire range of lifetimes is relatively 
long-lived on the time scales relevant to collider physics. 

Fortunately, some fraction of metastable charged or 
colored NLSPs, such as the squark, slepton, chargino, 
or gluino, will interact with and eventually stop within 
the material of the LHC detectors. A number of groups 
evaluated this stopping efficiency, as well as prospects 
for performing precision spectroscopic measurements on 
the ensuing late NLSP decays [l^. Recently a search 
for stopped gluinos performed by the CMS collaboration 
placed a bound of nig < 398 GeV for a stable gluino [31 ■ 
More generally, it is expected that CMS will effectively 
probe lifetimes of stopped particles in the range of 10^^ — 
10^ sec [13]. 

In the case of neutral NLSPs, such as the neutralino or 
sneutrino, stopping will not occur. That said, a sizeable 
fraction i/7CTNLSP of events may still decay within the 
length of the detector, L ~ 0(1 m), allowing for the 
possibility of lifetime probes in the range 10~^ — 10~^ 
sec at high luminosity. 

The key point is that there exists a precision correla- 
tion between the superpartner spectrum and the NLSP 
lifetime which, if verified, would provide very strong ev- 
idence for gravitino dark matter arising from freeze-in. 
Also, while Tr cannot be inferred accurately from col- 
lider measurements, precisely because freeze-in occurs on 
the near vertical part of the contours in Fig. [T] and Fig.[2l 
the reheating temperature will have a very strong upper 
bound at the order of magnitude level. 



III. CONCLUSIONS 

If the messengers of supersymmetry breaking are be- 
low the Planck scale, then the gravitino is the LSP and 
is thus a prime candidate for dark matter. We find that 
for a large range in Tr, gravitino dark matter is predom- 
inantly produced by freeze-in and is thus insensitive to 
Tfl. Furthermore, there is a direct correlation between 
the cosmological abundance of dark matter and the de- 
cay rate of the NLSP to gravitinos. The NLSP lifetime, 
given in Eq. ([8]) and shown as contours in Fig. [3l allows 
for a precision test at the LHC of the freeze-in origins of 
gravitino dark matter. Moreover, the reheat temperature 
can then be inferred, at least to within a couple of orders 
of magnitude, from Fig. [1] 

In [2^ we will provide additional motivations for grav- 
itino freeze-in through an investigation of the cosmology 
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of the QCD axino, the supersymmetric partner of the 
QCD axion. Our discussion will center on a "QCD ax- 
ino problem" which is entirely analogous to the gravitino 
problem but which occurs in a complementary region of 
Together, the combined axino and gravitino prob- 
lem completely exclude the possibility of a high Tr and 
thus much of parameter space in which gravitino pro- 
duction arises from thermal scattering. This indicates a 
robust tension between the axion solution to the strong 
CP problem and supersymmetry which strongly favors a 



low reheating temperature. 
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